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Abstract 

We present the applications of variational-wavelet ap- 
proach for computing multiresolution/multiscale represen- 
tation for solution of some approximations of Vlasov- 
Maxwell-Poisson equations. 

1 INTRODUCTION 

In this paper we consider the applications of a new nu- 
merical-analytical technique which is based on the meth- 
ods of local nonlinear harmonic analysis or wavelet analy- 
sis to the nonlinear beam/accelerator physics problems de- 
scribed by some forms of Vlasov-Maxwell-Poisson equa- 
tions. Such approach may be useful in all models in which 
it is possible and reasonable to reduce all complicated prob- 
lems related with statistical distributions to the problems 
described by systems of nonlinear ordinary /partial differen- 
tial equations with or without some (functional)constraints. 
Wavelet analysis is a relatively novel set of mathemati- 
cal methods, which gives us the possibility to work with 
well-localized bases in functional spaces and gives for the 
general type of operators (differential, integral, pseudodif- 
ferential) in such bases the maximum sparse forms. Our 
approach in this paper is based on the variational-wavelet 
approach from [1]-[10], which allows us to consider poly- 
nomial and rational type of nonlinearities. The solution has 
the following multiscale/multiresolution decomposition via 
nonlinear high-localized eigenmodes 
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which corresponds to the full multiresolution expansion in 
all time/space scales. 

Formula (jl|) gives us expansion into the slow part u'^ffj 
and fast oscillating parts for arbitrary N, M. So, we may 
move from coarse scales of resolution to the finest one 
for obtaining more detailed information about our dynam- 
ical process. The first term in the RHS of formulae (1) 
corresponds on the global level of function space decom- 
position to resolution space and the second one to detail 
space. In this way we give contribution to our full solu- 
tion from each scale of resolution or each time/space scale 
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or from each nonlinear eigenmode (Fig.l). The same is 
correct for the contribution to power spectral density (en- 
ergy spectrum): we can take into account contributions 
from each level/scale of resolution. Starting in part 2 from 
Vlasov-Maxwell-Poisson equations we consider in part 3 
the approach based on variational-wavelet formulation in 
the bases of compactly supported wavelets or nonlinear 
eigenmodes. 




Figure 1 ; Multiscale/eigenmode decomposition. 

2 VLASOV-MAXWELL-POISSON 
EQUATIONS 

Analysis based on the non-linear Vlasov-Maxwell-Poisson 
equations leds to more clear understanding of the collecti- 
ve effects and nonlinear beam dynamics of high inten- 
sity beam propagation in periodic-focusing and uniform- 
focusing transport systems. We consider the following 
form of equations ([11] for setup and designation): 
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The corresponding Hamiltonian for transverse single-par- 
ticle motion is given by 
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H{x,y,Px,Py,s) = -{p^ +py) + -[kx{s)x (5) 
+ky{s)y'^] + Hi{x,y,Px,Py,s) +ilj{x,y,s), 

where Hi is nonlinear (polynomial/rational) part of the full 
Hamiltonian. In case of Vlasov-Maxwell-Poisson system 



we may transform (2) into invariant form 

0. 
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3 VARIATIONAL MULTISCALE 
REPRESENTATION 

The first main part of our consideration is some variational 
approach, which reduces initial problem to the problem 
of solution of functional equations at the first stage and 
some algebraical problems at the second stage. Multires- 
olution expansion is the second main part of our construc- 
tion. Because affine group of translation and dilations is 
inside the approach, this method resembles the action of a 
microscope. We have contribution to final result from each 
scale of resolution from the whole infinite scale of increas- 
ing closed subspaces Vj: ■■■V-2 C V-i C Vb C C 
V2 C .... The solution is parameterized by solutions of 
two reduced algebraical problems, one is nonlinear and the 
second are some linear problems, which are obtained by the 
method of Connection Coefficients (CC)[12]. We use com- 
pactly supported wavelet basis. Let our wavelet expansion 
be 



t{x) + ^^Cjk-4^jk{x) (7) 
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If Cjk = for j > J, then f{x) has an alternative ex- 
pansion in terms of dilated scaling functions only f{x) — 
cj£(pji{x). This is a finite wavelet expansion, it can 

fez 

be written solely in terms of translated scaling functions. 
To solve our second associated linear problem we need 
to evaluate derivatives of f{x) in terms of (p{x). Let be 
(fi" = d"-(pi{x)/dx^\ We consider computation of the 
wavelet - Galerkin integrals. Let /"^{x) be d-derivative of 
function f{x), then we have /"^{x) = J^i'^ifei^)' ^^'^ 
values ip({x) can be expanded in terms of (p{x) 



(8) 



According to CC method [12] we use the next construc- 
tion. When N in scaling equation is a finite even positive 
integer the function (p{x) has compact support contained in 
[0, iV- 1]. For a fixed triple (di, dz, da) only some Ag^^^"^^ 
are nonzero: 2-N<£<N-2, 2 - N < m < 
N~2, \£-m\ < N-2. There are M = 3N^-9N + 7 
such pairs {£, m). Let A'^i*^'^^ be an M-vector, whose com- 
ponents are numbers A^^^^''^. Then we have the first re- 
duced algebraical system : A satisfy the system of equa- 
tions {d — di + d2 + ds) 

j^j^did2d3 _ 2i-dj^did2d3 (11} 
P 

By moment equations we have created a system of M + d+ 
1 equations in M unknowns. It has rank AI and we can ob- 
tain unique solution by combination of LU decomposition 
and QR algorithm. The second reduced algebraical system 
gives us the 2-term connection coefficients (d = di + ^2): 
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For nonquadratic case we have analogously additional lin- 
ear problems for objects (9). Solving these linear problems 
we obtain the coefficients of reduced nonlinear algebraical 
system and after that we obtain the coefficients of wavelet 
expansion (7). As a result we obtained the explicit time so- 
lution of our problem in the base of compactly supported 
wavelets. Also in our case we need to consider the exten- 
sion of this approach to the case of any type of variable 
coefficients (periodic, regular or singular). We can pro- 
duce such approach if we add in our construction additional 
refinement equation, which encoded all information about 
variable coefficients [13]. So, we need to compute only 
additional integrals of the form 



6,,(i)((^i)'*K2"i-fci)(<^2)''H2"t-fc2)dx, (13) 
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where are wavelet-Galerkin integrals. The coefficients 
A„i are 2-term connection coefficients. In general we need 
to find {di > 0) 
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For quadratic nonlinearities we need to evaluate two and 
three connection coefficients 
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where bij (t) are arbitrary functions of time and trial func- 
tions tfi , ip2 satisfy the refinement equations: 
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If we consider all computations in the class of compactly 
supported wavelets then only a finite number of coefficients 
do not vanish. To approximate the non-constant coeffi- 
cients, we need choose a different refinable function ips 
along with some local approximation scheme 

{Bif)ix) := FiA.f)M^'t - fc), (15) 

qGZ 

where F(,k are suitable functionals supported in a small 
neighborhood of 2^^k and then replace bij in (13) by 
Bebij{t). In particular case one can take a characteristic 



function and can thus approximate non-smooth coefficients 
locally. To guarantee sufficient accuracy of the resulting 
approximation to ( |T3[ ) it is important to have the flexibility 
of choosing (/^a different from (pi, <y92- In the case when D 
is some domain, we can write 



b^Jit) \d= Y1 Kit)XD{2't-k), 
0<fe<2* 
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where xd is characteristic function of D. So, if we take 
'P4 = Xd, which is again a refinable function, then the 
problem of computation of ( |l3| ) is reduced to the problem 
of calculation of integral 

V53(2'^t - k2)^p't{2''t - A:3)^2' - kA)dx (17) 

The key point is that these integrals also satisfy some sort 
of refinement equation [13]: 



2-\>'\H{k)^Y.b2k-iH{i), fi = di 
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This equation can be interpreted as the problem of comput- 
ing an eigenvector Thus, the problem of extension of the 
case of variable coefficients are reduced to the same stan- 
dard algebraical problem as in case of constant coefficients. 
So, the general scheme is the same one and we have only 
one more additional linear algebraic problem by which we 
can parameterize the solutions of corresponding problem in 
the same way. 

So, we use wavelet bases with their good space/time lo- 
calization properties to explore the dynamics of coherent 
structures in spatially-extended stochastic systems. After 
some ansatzes, reductions and constructions we give for 
(2)-(6) the following representation for solutions 

fc e 

where V^{s), Ui{z) are both wavelets or nonlinear high- 
localized eigenmodes and z = {x, y). 

Resulting multiresolution/multiscale representation for 
solutions of (2)-(6) in the high-locaUzed bases is demon- 
strated on Fig. 2. 
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Figure 2: The solution via multiscales. 
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